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Abstract 

We investigate the combined effects of electrostatic interactions 
and hydrodynamic interactions (HI) on the short-time rotational self- 
diffusion coefhcient in charge-stabilized suspensions. We calculate 
Dl as a function of volume fraction (f) for various effective particle 
charges and various amounts of added electrolyte. The influence of HI 
is taken into account by a series expansion of the two-body mobility 
tensors. At sufficiently small this is an excellent approximation due 
to the strong electrostatic repulsion. For larger 0, we also consider the 
leading hydrodynamic three-body contribution. 

Our calculations show that the influence of the HI on is less 
pronounced for charged particles than for uncharged ones. Salt-free 
suspensions are particularly weakly influenced by HI. For these strongly 
correlated systems we obtain the interesting result = Dq{1 ~ ar(t)^) 
for small 4>. Here DJj denotes the Stokesian rotational diffusion co- 
efficient, and ar is a positive parameter which is found to be nearly 
independent of the particle charge. The quadratic ^-dependence can 
be well explained in terms of an effective hard sphere model. 

Experimental verification of our theoretical results for is possible 
using depolarized dynamic light scattering from dispersions of optically 
anisotropic spherical particles. 



1 Introduction 

Until several years ago, the major interest in the dynamics of colloidal sus- 
pensions was focussed on the investigation of translational diffusion. The 
properties of translational collective and self-diffusion have been studied in 
great detail both experimentally Q and theoretically [^, ^, |^. Information 
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on the translational diffusion is contained in the dynamic structure factor 
S{q, t), which can be probed for a certain range of wavenumbers q and corre- 
lation times t using dynamic hght scattering from essentiahy monodisperse 
particles [P. 

More recently, considerable effort was made to investigate also the rota- 
tional diffusion in suspensions of spherically shaped colloidal particles, which 
possess an intrinsic optical anisotropy due to a partially crystalline struc- 
ture 1^, ^. To investigate the rotational particle diffusion, one has to resort 
to depolarized dynamic light scattering (DDLS), which is sensitive on the 
temporal self-correlations both in the particle positions and orientations. 
With this technique, the information on the short-time rotational diffusion 
is obtained from the autocorrelation function GvH{Q,t) of the horizontally 
polarized scattered electric field Q . From DDLS measurements of the first 
cumulant of GvHiQ,t), one can determine the short-time rotational self- 
diffusion coefficient D^. The diffusion coefficient depends on the direct 
potential interactions between the colloidal particles and on the indirect 
hydrodynamic interactions (HI) mediated by the suspending solvent. The 
latter interactions account for the fact that the velocity field, generated in 
the surrounding fluid by the motion of one particle, affects that of the other 
particles. 

By performing DDLS on index-matched suspensions of charged spher- 
ical particles made of a fluorinated polymer, Degiorgio et. al. were able 
to measure in great detail the concentration dependence of the rotational 
diffusion coefficient ^. They were mainly interested in suspensions with 
large salinity, where the particles essentially behave as hard spheres. Their 
experimental findings for have been compared with recent theoretical 
calculations by Jones designed for hard sphere suspensions. These 

calculations include also the approximative evaluation of the coefficient of 
the quadratic term in the virial expansion of in powers of the volume 
fraction (p. The agreement between theoretical and experimental results for 

is found to be satisfactory. Moreover, the small differences observed in 
the magnitude of the first and second virial coefficients are suggested to be 
partially due to residual electrostatic interactions, which are left over if the 
ionic strength of the added electrolyte is not large enough to completely 
screen the Coulomb repulsion 

In fact, the microstructure of a suspension of charged particles depends 
crucially on the amount of added electrolyte. For charge-stabilized sus- 
pensions, the general form of, e.g., the pair distribution function g{r) is 
qualitatively different from that of hard spheres at the same volume frac- 
tion [^. Due to the large spatial extend of the counterion cloud, the pair 
potential can have a range much larger than the physical particle diameter 
a. Therefore, appreciable pair correlations become already important at 
much lower values of the volume fraction (j) as in the case of hard spheres. 
For this reason, the concentration dependence of can be quite different 



2 



for uncharged and charged particles, particularly if the amount of added 
electrolyte becomes low. 

In this work, we focus on the combined effects of electrostatic and hydro- 
dynamic interactions on the short-time rotational self-diffusion in charge- 
stabilized suspensions of spherical particles. Starting from the generalized 
Smoluchowski equation for the positional and orientational degrees of free- 
dom, we calculate -D^ as a function of volume fraction (p for various effective 
particle charges and various amounts of added electrolyte. The hydrody- 
namic interactions are taken into account by a series expansion of the two- 
body mobility tensors in terms of the reciprocal interparticle distance, r~^, 
by including contributions up to order r~'^^. At sufficiently small cp, this in 
an excellent approximation due to the strong electrostatic repulsion, which 
renders configurations of nearly touching particles as extremely unlikely. 
For larger (j), we consider also the leading hydrodynamic three-body con- 
tribution Q, using Kirkwood's superposition approximation for the triplet 
correlation function, which is needed as the static input. 

We will show that the influence of the HI on is less strong for 
charge-stabilized particles than for uncharged ones. The influence of HI 
is particularly weak for suspensions in which all exess ions have been re- 
moved. For these strongly correlated systems, we find the remarkable result 
Dl = Dg(l— a.r0^) for small (j), where Dq is the rotational diffusion coefficient 
at infinite dilution. The parameter is found to be nearly independent of 
the particle charge for values large enough so that the hard core is com- 
pletely masked by the electrostatic repulsion. We can explain the quadratic 
^-dependence in terms of a simplified calculation based on an effective hard 
sphere model. 

The paper is organized as follows. In Sec. |2|, we briefly summarize 
salient relations of the theory of DDLS from optically anisotropic particles, 
which are relevant for obtaining the first cumulant of the depolarized field 
autocorrelation function. For later comparison, we recall in Sec. ^ the 
theory of short-time rotational diffusion and its main results with regard 
to suspensions of hard spheres. Sec. |^ gives the details of our calculations 
of D'g for charge-stabilized suspensions. It contains also a discussion of the 
range of validity of the various approximations used in this work. We present 
our results in Sec. ^ and discuss their meaning in terms of simple physical 
arguments based on an effective hard sphere model. Sec. ^ contains our 
final conclusion. 

2 Depolarized light scattering from anisotropic par- 
ticles 

In this section, we summarize some pertinent relations of the theory of depo- 
larized dynamic light scattering (DDLS) from optically anisotropic particles. 
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which are required for our further discussion. A thorough discussion of gen- 
eral DDLS properties has been given in [P, |6|, 0], which allows us to be rather 
brief in our explanations. 

Consider a suspension of identical spherical particles with cylindrically 
symmetrical optical anisotopy and in a situation in which the incident elec- 
tric field of the laser beam is linearly polarized perpendicular to the scatter- 
ing plane. The suspended particles are assumed to be nearly index-matched 
to the solvent so that the Rayleigh-Gans-Debye approximation is valid even 
for large volume fractions. The presence of optical anisotropy in the scat- 
terers gives then rise to a non-vanishing horizontally polarized part of the 
scattered electric field with magnitude 

EvHiq, t) = fiq)P i^-^j " (Y2,i{um + Y,^_,{uim e'^-'^^^'l (1) 

which is not affected by multiply scattered light. Here N is the number of 
particles in the scattering volume, q is the scattering vector with modulus 
q, and Ri{t) and ui{t) are the position vector and, respectivly, the unit 
vector in the direction of the optical axis of the Ith. particle. The l2,±i are 
the second-order spherical harmonics of index ±1, f{q) is the scattering 
amplitude of a particle, and /? = ay — a±, is the internal particle anisotropy, 
i.e. the difference in the polarizabilities parallel and perpendicular to the 
optical axis. 

In DDLS experiments, one measures the modulus of the temporal cor- 
relation function ^ 

GvHiq, t) = {EvHiq, 0)E*vh{<1, t)) (2) 

of the depolarized scattered field, where the bracket indicates for an ergodic 
system likewise a time or equilibrium ensemble average. 

On assuming that the orientation of a particle is decoupled from the 
particles translation, \GvH{Q,t)\ can be factorized as ^ 

\GvH{q,t)\ = Nf{q)^Fs{q,t)Fr{t), (3) 
15 

where 

Frit) = 4^(^2:1 (wi(o))y2,i(wiW)) (4) 

is the rotational self-correlation function and 

F,(g,t) = (e^'?-['«i(0)-«^Wl) (5) 

denotes the translational self-correlation function, also known in the litera- 
ture as the self-intermediate scattering function. The assumption of trans- 
lational and rotational decoupling which leads to Eq. (^) can be shown 
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in case of spherically symmetric interacting particles to be rigorously true 
at short times, i.e. to first order in t in a short-time expansion based on 
the generalized Smoluchowski equation |^. At longer times, the assumed 
decoupling is not strictly valid, but for hard sphere suspensions there is at 
least experimental evidence that deviations are small for all times ^. 
If HI are neglected, Fr{t) becomes exponential 

Frit) = (6) 

where Dq = ksT / {Sirrjo?) is the Stokesian rotational diffusion coefficient of 
a single particle of radius a, and r] is the shear viscosity of the suspending 
fluid. Eq. (^) arises from the fact that the orientational diffusion of particles 
with radial symmetric pair interactions are independent from each other if 
the HI are not considered. In the case of non-interacting spherical particles, 
Fs{q,t) also reduces to an exponential function , i.e. 

F,(g,t) = e-«'^o*, (7) 

where Dq = ksT / [Qn-qa) is the translational diffusion coefficient at infinite 
dilution. 

We focus now on the short-time behaviour of interacting colloidal par- 
ticles. The short-time rotational diffusion coefficient is defined by the 
short-time behaviour of Fr{t) as 

This measurable quantity contains the configuration-averaged effects of the 
HI on the short-time rotational diffusion of a spherical particle. As we will 
show in the next section, crucially depends an system parameters like 
the volume fraction the amount of added electrolyte and the effective 
particle charge. One finds -D^ — > -Dq in the limit — > 0. 

Notice that the short-time limit i ^ in Eq. (^) actually means that 
<C i <C with r]j = 0/(87rr7a^) being the relaxation time of the 

angular momentum of a colloidal sphere with moment of inertia 0. For uni- 
form spherical particles is ~ > where = m/ (Girrja) is the momentum 
relaxation time of a particle with mass m Q. It is ~ 10^*^5 for typical 
suspensions. Most dynamic light scattering experiments are restricted to 
correlation times t > 10~^s S> t|^. As a consequence, inertial effects arising 
from the linear and angular momentum relaxation of the particles are not 
resolved, so that in DDLS only the relaxation of the particle orientations 
and positions is probed. This fact allows for a coarse-grained description 
of Brownian motion on the basis of the generalized Smoluchowski equation, 
which describes essentialy this relaxation |3|, ^. On this level of description, 
HI can be considered to act instantaneously. 
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We further note that = r//3, where tj = (t'^/Dq is the socahed 

structural relaxation time, i.e. the time roughly needed for a non-negligible 
change of the direct interactions due to configurational relaxation. Typically 
one finds r/ ~ 10~^s such that the short-time regime <S t ^ r/ ~ I/-D0 
is well seperated from the long-time regime t ^ tj. For the latter regime, 
memory effects are important and lead to deviations in \GvHiQ,t)\ from 
purely exponential decay. 

3 Short— time rotational diffusion 

For calculating D^, it is necessary at first to specify the time evolution of 
the suspended particles. Since we restrict ourselves to correlation times 
i S> r]^, the suspension dynamics is properly described by the generalized 
Smoluchowski equation. In Smoluchowski dynamics, Fr{t) can be expressed 
as 

Frit) = 47r(y2:i(«i(0))e^^*y2,iK(0))), (9) 

where 

= ^ { (v. - ■ (d% ■ V, + D% ■ L,) + 

ii-(£>l|-V,+£>:}'-i,)|, (10) 

is the adjoint (or backward) Smoluchowski operator [^, ^, Here Lj = 
Uj X -^j- is the gradient operator in the compact space of orientations of the 
jth particle, with Uj being the unit vector pointing in the direction of its 
optical axis. In writing Eq. (Il0| ) it has been assumed that the total energy 
$ = ^{R^) of the N particles depends only on the center-of-mass positions 
R^ = [Ri, ■ ■ ■ Rn) and not on the particle orientations. We further assume 
<^ to be pairwise additive, i.e. 

1 ^ 

'^{R'') = -Y.'u{Rm), (11) 

k,l=l 

where u{r) is a spherically symmetric pair potential, and Rm = \Rki\ with 
Rki = Rk — Ri- The prime indicates that the term k = I \s excluded from 
the sum. The diffusivity tensors -D^j, with a,b £ i}, embody the HI 
between the particles by coupling the forces = {Fi, . . . Fj^j) and torques 
j^N _ ^ji^^ ^ ^ ^ Tjv) acting on the particles to their translational velocities 
= (Vi, . . . Vn) and angular velocities u'^ = (ui, . . .u^): 

_ ^/£>"(i?^)£>*^(i?^)\ (F^ 
\ujN ) ~ keT \D'-^R^)D'-'-{R^) J 
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For spherical particles, the diffusivity tensors depend only on the position 



vectors. The tensors = Dfj /{ksT) are called mobility tensors. 
Using Eqs. (^) and (^), one can express into the form 



d: 



and the dimensionless diffusion coefficient is given by 



h: 



1 



{T^DYiiR^)). 



3DI 



(13) 



(14) 



Here TrD\\ denotes the sum on the diagonal elements of the tensor DW- 
Due to the many-body nature of the HI, it is not possible to perform an 
exact evaluation of H^, valid at all particle concentrations. For small vol- 
ume fractions however, when the mean interparticle distance gets sufficiently 
large, it becomes possible to obtain a good approximation for by consid- 
ering only two-body and, to leading order, three-body contributions to the 
HI. 

By using a rooted cluster expansion, Degiorgio et.al. have shown that 
the normalized short-time rotational diffusion coefficient can be expressed 
as a series S 



(15) 



The coefficient H^i of the linear term is expressable in term of integrals 
involving only hydrodynamic two-body contributions and the radial distri- 
bution function g{r). Explicitly, Hg, is given by |^ |9| 



H 



si 



(16) 



and this expression involves two scalar mobility functions ctyKr) and PlKr) 
depending on the interparticle distance r. These functions are calculable by 
means of a series expansion in even powers of (a/r) [|l^]. We only quote the 
leading terms 




15 /a 



39 /a 



+ 




Inserting Eqs. <^^^) in Eq. (jl^) leads to 



H 



si 



dt g{t) 



2 2 V 



(17) 
(18) 

(19) 



with t = r/a. 
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The second coefficient in the expansion of Eq. ( [TBD is for more dif- 
ficult to evaluate since it involves hydrodynamic three-body contributions. 
Using the method of reflections, Jones (in Ref. [^) was able to calculate the 
leading term in the far-field expansion of the irreducible three-body mobil- 
ity tensor. Considering only this term, is approximated by a three-fold 
integral 

225 /"^ /"^ 

^l2 = -^ X Jo 7_/^l5^'^(*12,il3,6)/(il2,tl3,6), (20) 

which involves the static three-particle distribution function g^^^ in depen- 
dence of ti2 = 2a/Ri2, ti3 = 2a/i?i3, and ^ = R12 ■ Ris/iRuRn). The 
somewhat lengthy expression for the function /(ti2iii3iCi) is given in Eqs. 
(36-37) of Ref. and will not be repeated here. 

For dilute to moderately concentrated suspensions, it is sufficient to con- 
sider only the leading three terms in the expansion Eq. (15), with H^^ and 



HI2 expressed through Eqs. (|16D and (|20|). It should be stressed that Eq. 
(15) is, up to this point, not just a virial expansion of -ffj in powers of 



since g{r) and g^^'^{r,r^) are also (^-dependend. To proceed further, one 
needs to specify the static distribution functions, which themselves depend 
crucially on the form of the pair potential u{r). It is at this point where the 
important differences in the short-time behaviour between charge-stabilized 
suspensions and suspensions of hard spheres come from. 

For later comparison consider first a dilute suspension of colloidal hard 
spheres (typically with cp < 0.1). It is then allowed to use in Eqs. (|16|) and 
(|20|) for consistency the leading order virial expansion of the static distribu- 



(3) 

tion functions gnsi''') ^-nd 5)75 (^;^') to first and zeroth order, respectivly. 
If the exact numerical input for the two-body functions aii{r) and /^[[(r) 
is employed, the following truncated virial expansion for is obtained |^ 

= l- 0.630</.- 0.6702 + 0(</,3). ^21) 

The radial distribution function gnsi''') has its maximum at contact, i.e. 
9Hsi^ = <7^) > 1, and the contact value is a monotoneously increasing 
function of the volume fraction in the liquid regime 4> < 0.49. In the case of 
hard spheres, it is therefore necessary to include many terms in the expansion 
of the mobility tensors in powers of a/r. In fact, if the two-body scalar 
functions in Eq. ([l^ ) are approximated by series expansions up to 0{r~^^), 
we would get a value of —0.578 for the first virial coefficient, significantly 
different from the exact numerical result —0.630. 



4 Calculation of for charge— stabilized particles 

In the present work, we study the short-time rotational diffusion of suspen- 
sions of charge-stabilized particles. It is crucial to note in contrast to hard 
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sphere suspensions that charge-stabihzed systems at sufficiently low ionic 
strength exhibit even at small volume fractions spatial correlations with pro- 
nounced oscillations in g{r) (cf. Fig. For this reason, it is in general 
not possible for charge-stabilized suspensions to use a virial expansion for 
the static distribution functions. Consider for example the zeroth density 
limit goir) = exp [— /3u(r)] of the radial distribution function. Due to the 
long-range nature of the electrostatic pair forces, extremely small values 
of (j) are needed for g(r) to approach this limit. As a consequence, g{r) 
has to be calculated from appropriate integral equation schemes or from 
time-consuming computer simulations. In other words, charge-stabilized 
suspensions at small to moderately large (j) can be regarded as dilute with 
respect to the suspension hydrodynamics (which allows to use a rooted clus- 
ter expansion truncated after the three-body contributions), but they are 
"concentrated" as far as the microstructure is concerned. 

In performing the numerical integration in Eq. (0), we have decided 
to calculate g{r) using the rescaled mean spherical approximation (RMS A) 



[11 1, mainly for its numerical simplicity but also since it has been found to 



be an efficient fitting device of experimentally determined structure factors 

i- 

Our RMSA results for g{r) are based on the effective macroion fluid 
model of charge-stabilized suspensions. In this model, the effective pair 
potential u{r) between two charged particles is described as a sum of a 
hard-sphere potential with diameter a and a screened Coulomb potential 

Here /3 = l/{kBT)^ Lb = /3e^/e is the Bjerrum length, e the elementary 
charge, e the dielectric constant of the suspending fluid, and Z the effective 
charge (in units of e) of a colloidal particle. The equation 

=47rLB[n|Z| + 2n,] (23) 

defines the Debye-Hiickel screening parameter k, where n and Ug are the 
number densities of colloidal particles and of an added 1-1 electrolyte, re- 
spectivly. It is assumed that the counterions are monovalent. 

The leading-term expression Eq. (|2^) for has the triplet distribution 
function g^^\r,r^) as input. There are, to the best of our knowledge, so far 
no manageable numerical schemes available which provide decent approx- 
imations for g^^\r,r^) in the effective macroion fiuid model. To calculate 
the three-body contribution given by Eq. ( ]20|), we therefore use for 
simplicity Kirkwood's superposition approximation [1^ 

g^'Hr,r')=g{\r-r'\)g{r)g{r'), (24) 
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with g{r) calculated again with the RMSA scheme. The three-fold integra- 
tion is performed using a Monte Carlo method. 

We have stressed that in the case of hard spheres many terms need to 
be considered in the far-field expansions of the hydrodynamic mobility ten- 
sors, since configurations of nearly touching hard spheres are very likely. The 
situation is quite different for charge-stabilized suspensions: the strong re- 
pulsion of the electric double layers keeps the particles apart from each other, 
resulting in a very small probability of two or more spheres getting close to 
each other. Indeed, the g{r) for charged particles remains essentially zero 
for particle seperations comparable to the Debye-Hiickel screening length 
(cf. Fig. I). 

Therefore, in charge-stabilized suspensions at sufficiently low ionic strength, 
it is possible to account only for the first few terms in the far-field expan- 
sions of the mobility tensors. In our calculations of -ffji, we include terms up 
to O ((a/r)20) in the expressions of the two-body mobility functions |lO| ]. 
Our results will show that only the first three terms up to O (^r~^^) are 
significantly contributing to D^. 

5 Results and discussion 

In this section we present our numerical results for the normalized short- 
time rotational diffusion coefficient of charge-stabilized suspensions as 
a function of volume fraction (p, effective charge Z, and concentration of 
added 1-1 electrolyte rig. The system parameters used in our calculations 
are typical for partially crystalline spherical particles made of a fluorinated 
polymer and dispersed in an index-matching solvent mixture of water with 
20% urea: e = 87.0, T = 294K, and particle diameter a = 900A [|13|. The 
effective charge chosen in our RMSA calculations of g{r) is Z = 500 for 
those results where Z is kept constant. If not stressed differently, we have 
included in our calculations of two-body contributions up to order r"^" 
together with the leading three-body contribution. 

In the following, we will show that the volume fraction dependence of 
Hg for charged particles is qualitativly different from that observed for hard 
spheres. The differences are most pronounced for deionized charge-stabilized 
suspensions, where essentially all excess electrolyte has been removed. Thus, 
we first concentrate on suspensions with Ug = 0. 

5.1 Deionized suspensions 

Fig. m displays our results for iJJ = D^/Dq as a function of (crosses). For 
comparison, the corresponding result, Eq. for hard sphere suspensions 

is also shown in the figure. Evidently, the influence of the HI on HI is 
less pronounced for charged particles than for uncharged ones. We will 
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later show that salt-free (i.e. deionized) suspensions are particularly weakly 
influenced by HI. 

A best fit of our results for HI in deionized suspensions gives the inter- 
esting result 

i/J = l-a^0^ ~ 1.15, (25) 
i.e. a quadratic (/)-dependence! This should be contrasted with the corre- 



sponding expression Eq. (21) for hard spheres, where the linear term in 



gives the dominant contribution to if (/) < 0.15. The quadratic 



dependence in Eq. (25) is valid up to surprisingly large volume fractions, 
typically up to ~ 0.25. Moreover, the coefficient ar is found to be nearly 
independent of the particle charge, for values of Z large enough so that the 
hard core of the particles remains completely masked by their electrostatic 
repulsion (i.e. typically for Z > 200). This fact is illustrated in Fig. ^, 
which shows results for Hg{(j)) for various values of Z. All graphs in this 
figure can be fitted by the functional form Eq. ([25|), with numerical values 
of ar close to 1.15. 

A physical explanation for the, as compared to hard sphere suspensions, 
weaker influence of the HI on in case of charged particles can be obtained 
from Fig. [l|. In this flgure, RMSA results for g{r) are shown for different 
values of (j). These graphs for g{r) illustrate the pronounced interparticle 
corrrelations prevailing in deionized suspensions down to very small volume 
fractions. Notice also from this flgure that the strong electrostatic repulsion 
gives rise to a "correlation hole" centered around each particle, i.e. a spher- 
ical region usually extending over several particle diameters with (nearly) 
zero probability for finding another particle. The size of this correlation hole 
increases with decreasing </>. As a result, there is at small (j) only a weak hy- 
drodynamic coupling between the rotational motion of two or more spheres, 
and the deviations of from its value -Dq infinite dilution becomes quite 
small. To quantify this point, we have plotted in Fig. || the two-body scalar 
mobility function from Eq. (|^), —S-rrrja^ [oiiiii^) + 2/9ii(r)] , versus the re- 
duced interparticle distance. According to Eqs. (|T7|-p^) and Fig. ||, this 
function is a rapidly decaying function of r. Obviously, the mobility function 
is very small at those values of r where g{r) is different from zero, with the 
consequence that the value of the integral in Eq. ( |l6|) is small indeed. 

We present now an intuitive physical explanation for the quadratic (^-de- 
pendence of Dl, based on a characteristic scaling property of the principal 
peak of g{r), combined with a crude approximation for g(r) which incor- 
porates this scaling property. In this approximation, the realistic g{r) is 
replaced by a step function 

gEHsir) = e(r - 2aEHs), (26) 

where uehs > a is an effective hard sphere radius, which accounts grosso 
modo for the electrostatic repulsion between the particles. We refer to this 
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simplified description of the real pair structure as the effective hard sphere 
(EHS) model. Using this model for g{r) together with the far-field expan- 
sions Eqs. (17-15) of the two-body mobility functions, it is easy to calculate 



the leading terms in Eq. (^). The result is an expansion 

Hl^ = -— - - ^x-^ + 0(x-9) (27) 

'^i 16 64 64 ^ ^ ^ ' 

in terms of the ratio of the effective radius to the actual particle radius, 
X = aEHs/cL > 1. At small (/>, one can neglect the three-body contribution 
-f^s2<^^ as compared to Hlicf), and is then approximated in the EHS model 
as 



h: = 1 



16 64 64 ^ ^ ^ 



(28) 



In the next step, we need to specify the effective radius. There are various 
ways to estimate the magnitude of the effective radius in terms of the param- 
eters of the actual system A reasonable choice adopted here and shown 
to be very useful in earlier applications of the EHS model, is Idj^ug — '^m* 
Here rm is the position of the principal peak of the actual g{r), as calculated 
using the RMSA. A typical RMSA-g'(r) and the corresponding EHS-g(r) 
are depicted in Fig. 

For deionized suspensions with completely masked hard-core repulsion, 
it is well known that and hence x obeys the following scaling property 
in terms of (p 

x (X rm (X <p 3 . (29) 



This behaviour is due to the strong electrostatic repulsion, which leads to 
r^n ^ f = a(47r/3(/))~^/^, where f denotes the geometrical average distance 
between neighboring particles. 

In very diluted suspensions, x becomes so large that it is sufficient to 
take into account only the leading term in Eq. (|2^), which arises from the 
leading term proportional to in the expression Eq. (|18|) of /?[[(r). The 
EHS model predicts for this case a quadratic (?f)-dependence, since from Eqs. 



(ID and (|2|) it follows 

5 



h: = 1 



-x~^ + 0{x~'] 



1- A</.2 + O(0*'/^), (30) 



with A determined to A ~ 0.60 if is approximated by f. Hence we can 
conclude that the simple EHS model is for small (j) in qualitative agreement 
with our more refined numerical results for based on a more realistic 
g{r). The large deviation between A and the coefficient in Eq. ( p5|) arises 
from the "fine structure" in g{r), not captured by the EHS model. 

Since is determined by the volume fraction only, the EHS model 
further predicts to be independent of the particle charge, provided that 
Z is large enough (i.e. > 200) for rm oc to be valid. This prediction 
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of the EHS model is again in good agreement with our numerical finding 
for the ^-dependence of the coefficient a^. in Eq. (|25|). That the two-body 
contribution to i/J, i.e. -ffji, is indeed not sensitive to changes in Z follows 
from the observation that the principal peak of g{r) becomes higher and 
narrower on increasing Z at constant (f), but its position is nearly constant. 
The overlap region between g{r) and + '^Prii''")] which determines 

remains therefore also nearly constant. 

The EHS model suggests that the quadratic (/>-dependence of HI in 
our numerical results arises for small (j) from the leading term in the far- 
field expansion of the two-body mobility function [aij^r) + 2/3[j"(r)]. For a 
numerical check of this assertion consider Fig. |5|, which shows our result for 
HI if only the lowest order two-body contribution of 0{r~^) is considered 
(dashed line) . In comparison we show the result of the full calculation (full 
line), which was already displayed in Fig. ^, and which accounts for the 
leading three-body term and all two-body contributions up to 0(r~^°). 
We notice that both lines nearly superimpose on each other even at larger 
volume fractions (j) > 0.1, where the three-body term and the higher order 
two-body terms are expected and found to give non-negligible contributions 
to -D^. However, at larger (j) there is a fortuitous cancellation between the 
leading three-body contribution Hg2'p'^ and the higher order (i.e. 0{r~^)) 
two-body contributions to H^icp, which leaves the lowest order two-body 
contribution as the most significant term even at larger (p. 

Fig. H includes as the dotted line the result obtained for H^ if the three- 
body contribution is neglected but all two-body contributions up to 0{r~'^^) 
are accounted for. We obtain the same result if only two-body contributions 
up to 0{r~^^) are considered. Even the difference observed in the results 
for HI including two-body terms up to 0{r~^) and 0{r~'^^), respectively, is 
very small (cf. Fig. |5|). As a conclusion, in charge-stabilized suspension it 
is justified to use a truncated far-field expansion of the mobility functions 
even at larger (j). 

The fortuitous cancellation between the leading three-body term and 
the higher order two-body contributions observed in Fig. |5| could have 
been anticipated from the EHS model. Using the first three terms in the 
expansion of Eq. (|2^), we obtain the following expansion in cf) by using Eq. 



16 64 64 ^ ' 



-A(f^ - B(pl -C(l)^ +o((t)^Y (31) 



with constants A,B,C > 0. For an estimate of B and C, let us approximate 
rm again by f. This leads to i? ~ 0.41 and C ~ 0.21. We do not need 
to consider higher order terms in Eq. (|3l[), since these do not contribute 
significantly to i/Jj^. 
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Eq. ( pO| ) for yields together with Kirkwood's superposition approx- 
imation, QEHsir) (cf. Eq. 



s2' 



|26|)), and Eq. ( |29D the following result 
= 0.339 = D(j)^ 



(32) 



for the three-body contribution to H^. This result includes our finding 



(33) 



which states that the three-body contribution {H12)ehs proportional to 
(j). The numerical value {H12)hs ~ 0-339 was first obtained by Degiorgio 
et. al. II using a Monto Carlo integration method. Assuming r^n — r, 
D is determined as ~ 0.65. It is now evident that there is a partial 
concellation between the negative contribution —B(j)^^^ — C(j)^^^^ and the 
positive contribution Dcj)^, leaving the term —AcjP' as the most significant 
contribution to H^. 



5.2 Suspensions with added electrolyte 

We discuss now the dependence of D^/Dq on the amount od added elec- 
trolyte. In this context, it should be noted that most "hard-sphere" sus- 
pensions studied so far by DDLS are in fact suspensions of charged particles 
with a large amount of salt added to screen the electrostatic repulsion (cf. 
e.g. 1^, ^). Our results for Hl{(p) for various amounts of added 1-1 elec- 
trolyte ranging from = to = 13mM are displayed in Fig. |^. It is 
noted from this figure that addition of electrolyte leads to a decrease of H^, 
i.e. the rotational diffusion of the particles becomes more affected by the 
HI with increasing ionic strength. This is due to the fact that the system 
gradually transforms with increasing into a hard-sphere-like dispersion. 
The behaviour of can be qualitatively explained on the basis of the EHS 
model by noting for increasing that the effective radius uehs = de- 
creases and eventually approaches the physical diameter a for very large n<j. 
It follows from Eq. ( p^ ) that becomes smaller with decreasing x oc uehs- 
Going beyond the EHS model, we can learn more about the ionic strength 
dependence of by considering Fig. which shows RMSA results for 
g{r) for a suspension of volume fraction (j) = 0.01 and various amounts of 
added electrolyte. The g{r) of the deionized suspension has pronounced os- 
cillations, indicating strong particle correlations. These oscillations become 
damped out as electrolyte is added, and the peak position is shifted 
towards smaller values. For the system under consideration, very large 
amounts of electrolyte (i.e. ^ 13mM) are needed to reach the hard- 
sphere limit, where the RMSA-g'(r) becomes equal to the corresponding 
Percus-Yevick solution |jl2|. Consider, e.g., the case = ISmM in Fig. 
which corresponds to na ~ 33 and f3u{r = a^) ~ 6.0. The value of the pair 
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potential at contact distance is still large enough to create a small but non- 
negligible correlation hole with zero probability of finding another particle. 
In fact, it is well known |8|, ^ that the value of H^, and hence the value of 
Dl, is extremely sensitive to the behaviour of the radial distribution func- 
tion near touching. This is clearly seen from the fast decay of the two-body 
scalar mobility function — Svrr^a^ [aii{r) + 2/?[j"(r)] shown in Fig. 0. 

The scaling property oc is invalid for larger amounts of elec- 
trolyte, when the particle diameter becomes a second physically relevant 
length scale besides the mean particle distance f. The behaviour of Hl((p) 
at small (j) changes then gradually with increasing Ug from a quadratic to 
a linear (/)-dependence. At very large ^ WmM, g{r) exhibits only tiny 
changes in its form when (j) is increased. Hence gets more and more 
independent of 0, and Hl[(j)) = 1 + Hl^{(f))(j) + 0{(tP') becomes linear in (j). 
A deficiency of our calculations is the fact that we do not obtain the exact 
hard-sphere result {Hl^^jg = —0.630 in the limit — > oo. Instead the 
limiting value {Hl^)^g = —0.578 is reached, since we use a series expansion 
up to terms of order r^"^^ in calculating -ffJi. However, this truncated ex- 
pansion is sufficiently good for charge-stabilized suspensions with moderate 
amounts of added electrolyte, where x > 1.5 holds (cf. js], [l0t|). 

6 Conclusions 

In the present work, we have investigated the combined effects of electro- 
static interactions and hydrodynamic interactions on the short-time ro- 
tational diffusion coefficient in monodisperse suspensions of charge- 
stabilized colloidal particles. 

On the basis of the one-component macrofluid model, which neglects 
electroviscous effects arising from the dynamic distortion of the electric dou- 
ble layer around a particle, we have calculated iJJ = D^/Dq for various 
systems as a function of volume fraction, effective particle charge, and ionic 
strength. For these calculations we have used a rooted cluster expansion 
derived by Degiorgio et. al. Q, combined with far-field expansions of the 
two-body and three-body mobility functions. We have given evidence that 
this is a good and well founded approximation in case of charge-stabilized 
suspensions, as long as the amount of added electrolyte is not very large. 

The short-time rotational diffusion of charged particles is less affected 
by the HI than the rotational diffusion of the uncharged ones. For deionized 
suspensions of highly charged particles we have found a quadratic volume 
fraction dependence of the form = 1 — arcf^, with a coefficient a,. ~ 
1.15, which is nearly system independent. The quadratic (^-dependence 
extends up to rather large volume fractions due to the fact that the leading 
three-body term is essentially cancelled by the two-body contributions of 
0{r~^). We have further shown that the qualitative behaviour of can 
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be understood in the framework of a simplified EHS model. 

Incidently, dilute deionized suspensions are peculiar also with respect to 
the 0-dependence of the translational short-time self-diffusion coefficient 
D\. For this quantity, one obtains a non-analytic concentration dependence 
of the form D\/Dq = 1 — atcf)^/'^, with at — 2.5 The fractal exponent 
4/3 follows also from the EHS model applied to translational self-diffusion. 

Upon addition of salt, the quadratic (/)-dependence of changes grad- 
ually to a (/)-dependence typical for hard spheres, where the term linear in (f) 
dominates for (p < 0.1. Very large amounts of salt are required for the sys- 
tems investigated in this work to reach the hard-sphere limit. This finding is 
due to the fast decay of the two-body mobility functions with increasing r, 
which renders the two-body term Hgi extremly sensitive to distances close 
to the contact distance of two particles. 

Degiorgio et. al. Q report a small increase of the measured Dl of hard 
spheres at the freezing transition where </>/ = 0.49. The theoretical argument 
to explain this enhancement is the slightly increased free volume per particle 
in the crystal phase. The resulting slightly increased interparticle distance 
turns the HI somewhat less important. We expect that a similar effect 
should occur for the freezing transition in charge-stabilized suspensions, 
which occurs at much smaller (p. In our calculations we have used pair 
distribution functions g{r) calculated with a method (RMS A) designed only 
for the fluid phase. Therefore, we can not account for the effect of the 
freezing transition on the short-time rotational diffusion. Nevertheless, we 
anticipate that our calculations give reasonable results even for the colloidal 
crystal phase, since our arguments concerning the scaling property of the 
peak position of g{r) should also apply to the crystal phase. Beside that 
we have shown that Hg is not very sensitive on the details of the shape 
of g{r) near r^. For these reasons, we expect only a small decrease of 
the coefficient at the freezing transition. Therefore we have shown in 
this paper calculations of up to moderatly large (/>, although we know 
from the Hansen- Verlet criterion that the systems studied here should be 
crystalline already at small (p. 

To our knowledge no experimental or computer simulation data of 
for deionized charge-stabilized suspensions have been published so far. Be- 
cause of the interesting differences in the rotational diffusion between sus- 
pensions of charged and uncharged particles, it would be worth to check the 
general predictions of our calculations against experiments and computer 
simulations. We finally mention that recent DDLS measurements of in 
deionized suspensions of charged fluorinated polymer particles, initiated by 
our work, compare favorable with our calculations ||l3| . 

References 



16 



[1] p. N. Pusey, In "Liquids, Freezing and Glass Transition: 11", pp. 763- 
942, Les Houchcs Sessions 1989, edited by J. -P. Hansen, D. Levesque, 
and J. Zinn- Justin, editors. North Holland, Amsterdam, 1991. 

[2] W. Hess and R. Klein, Adv. Phys., 32, 173 (1983). 

[3] R. B. Jones and P. N. Pusey, Annu. Rev. Phys. Chem., 42, 137 (1991). 

[4] G. Nagele, Phys. Rep., in press (1996). 

[5] V. Degiorgio, R. Piazza, and T. Bellini, Adv. Coll. Int. Sci., 48, 61 

(1994) . 

[6] V. Degiorgio, R. Piazza, and R. B. Jones, Phys. Rev. E, 52, 2707 (1995). 

[7] B. J. Berne and R. Pecora, Dynamic Light Scattering, John Wiley &; 
Sons, New York, 1976. 

[8] R. B. Jones, Physica, A 150, 339 (1988). 

[9] B. Cichocki and B. U. Felderhof, J. Chem. Phys., 89, 1049 (1988). 

[10] R. B. Jones and R. Schmitz, Physica, A 149, 373 (1988). 

[11] G. Nagele, M. Medina-Noyola, R. Klein, and J. L. Arauz-Lara, Physica, 
A 149, 123 (1988). 

[12] J.-P. Hansen and I. R. McDonald, Theory of Simple Liquids, 2nd edi- 
tion, Academic Press, London, 1986. 

[13] F. Bitzer, T. Palberg, and P. Leiderer, University of Konstanz, private 
communication . 

[14] G. Nagele, B. Mandl, and R. Klein, Progr. Colloid Polym. Sci, 98, 117 

(1995) . 



17 



Figures 



2.0 



(|)=0.001 ^, 



--■ (|)=0.003 ■' 



-- (|)=0.01 



g(r) 



1.0 



0.0 



1.0 



3.0 



5.0 



r/o 



7.0 



9.0 



Figure 1: RMSA results for the radial distribution function g(r), modelling 
deionized suspensions at various volume fractions as indicated in the 
figure. System parameters are: a = 900A, Z = 500, e = 87.0, T = 294K, 
and Ug = 0. The full line is the graph of the two-body scalar mobility 
function —Sirrja^ [q^ii(^) + ^/^[[(r)], including terms up to 0{r^'^^) (cf. Eq. 

dD). 




Figure 2: Results for the normalized short-time rotational self-diffusion 
coefficient = D^/Dq for deionized systems, in comparison with the cor- 
responding hard-sphere result of Degiorgio et. al. [6]. Best fit has quadratic 
0-dependence: = 1 — 1.15(/>'^. System parameters as in Fig. |^. 
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Figure 3: Normalized short-time rotational self-diffusion coefficient HI for 
various values of the effective charge Z. All system parameters besides Z as 
in Fig. ||. 



g(r) 




Figure 4: Typical RMSA-g(r) (dashed line), and corresponding radial 
distribution function QEHsir) of the effective hard-sphere (EHS) model with 
'^ciEHS = fm- Note the scaling property oc <^~^/^, valid for deionized 
suspensions. 
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Figure 5: Dependence of HI on various two-body contributions, and on 
the leading three-body contribution. The various cases are indicated in the 
figure. Parameters as in Fig. ||. 




Figure 6: Volume fraction dependence of iJJ for various amounts of added 
1-1 electrolyte, as indicated in the figure. System parameters as in Fig. ^. 
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Figure 7: RMSA results for g(r) for (j) = 0.01, and for various amounts of 
added 1-1 electrolyte. Other parameters as in Fig. |l[. Full line: two-body 
mobility function — Svrr/a^ [Q^n('') + 2/3[[(r)] as in Fig. |^. 
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